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Summary 


Results are presented for an experimental study of the response of inertial and optical 
wind-tunnel model attitude measurement systems in a wind-off simulated dynamic 
environment. This study is part of an ongoing activity at the NASA Langley Research Center 
to develop high accuracy, advanced model attitude measurement systems that can be used 
in a dynamic wind-tunnel environment. This activity was prompted by the inertial model 
attitude sensor response observed during high levels of model vibration. The inertial sensors 
cannot distinguish between the gravitational acceleration and centrifugal accelerations 
associated with wind-tunnel model system vibration which results in a model attitude 
measurement bias error. Significant bias errors in model attitude measurement were found 
for the measurement using the inertial device during wind-off dynamic testing of a model 
system. The amount of bias present during wind-tunnel tests will depend on the amplitudes 
of the model dynamic response and the modal characteristics of the model system. 
Correction models were developed and used to predict the vibration-induced bias errors for 
the inertial device. A correction model based on the measured vibration amplitudes and 
modal characteristics of the model system predicts the bias error to a high degree of 
accuracy for the vibration modes characterized in the simulated dynamic environment. The 
optical system results were uncorrupted by model vibration in the laboratory setup. The 
optical model attitude measurement system and correction methods for the inertial model 
attitude measurement system are to be validated during actual wind-tunnel tests. 

Introduction 

The predominant instrumentation used to measure model attitude or angle of attack (AOA) in 
wind-tunnel testing at NASA Langley Research Center (LaRC) is the servo-accelerometer 
device described in reference 1 . For quasi-static conditions, this inertial sensor provides a 
model attitude measurement with respect to the local gravity field to an accuracy of ±0.01 ° 
over a range of ±20°. However, during wind-tunnel testing, dynamic pitching and yawing 
motion of the model creates centrifugal acceleration which is sensed by the onboard inertial 
AOA device and results in a bias error in the model attitude measurement. The amount of 
bias present during wind-tunnel tests will depend on the amplitudes of the model dynamic 
response and the modal characteristics of the model system. The problem of the inertial 
device sensitivity to model vibrations is briefly discussed in reference 1 . 

The National Transonic Facility (NTF) [2] is a transonic wind-tunnel located at NASA LaRC 
which has the capability fortesting models at Reynolds number up to 140 million at Mach 1 
and dynamic pressure up to 7000 pounds per square foot. Severe model vibrations have 
been encountered on a number of models since the tunnel began operation. References 3 
through 6 document studies of model and model support vibrations in the facility. Reference 
7 documents results of an experimental study conducted in 1 993 at the NTF to study the 
inertial AOA sensor response to a simulated dynamic environment. The experimental study 
[7] clearly established that AOA bias error is due to centrifugal forces associated with model 
vibration. During wind-off dynamic tests, bias errors over an order of magnitude greater than 
the desired device accuracy of 0.01 degree were measured. The bias error was found to be 
dependent on the vibration mode and amplitude. The study revealed the complexity of the 



problem when multiple vibration modes were present involving both pitch and yaw motions 

— er f r eC, ' 0n T 18 ' was devised from the P roblem Physics which gave good 
estimates of bias error when compared with test results. y 
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Inertial sensors have other problems besides the possible bias errors caused by model 

naUfrp TheS ® Pr ° b *'T S include: (1) retirements for special handling due to their fragile 
anH S 2 requ,r ® ments for multiple accelerometers for the measurement of pitch and rol? 
d (3) space and wiring requirements in the model fuselage. For these reasons (in addition 
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other tunnels) have been sought. Much of the interest has been in optical techniques which 
are generally assumed to be immune from the dynamic problems inherent in inertial devices 
However, even though an optical technique may not be theoretically corrupted by model 
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foings as variations in irradiance, camera/laser mounting system vibrations, and aliasing 
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Symbols List 


AOA 

angle of attack 

A?/ 

filtered AOA signal from inertial device 


peak acceleration for /#* vibration mode 

a n (t) 

time dependent normal acceleration 

a t (t) 

time dependent tangential acceleration 

a*(0 

time dependent longitudinal acceleration 

Aunf (0 

time dependent unfiltered AOA signal 

9 

gravitational constant 

Hz 

Hertz 

m 

number of included modes 
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q(t) time dependent inertial AOA output 

r current mode number 

Sqq power spectrum of q(t) 

Syy power spectrum of y(t) 

Sqy cross spectrum of q(t) and y(t) 

t time in seconds 

v r (f) time dependent velocity for r*h mode 

V r peak velocity for r^ 7 mode 

x(t) time dependent lateral motion (via balance or accelerometer output) 

y(t) integral-square of x(t) 

y r (t ) time dependent displacement for ft* mode 

Y r peak displacement for r*h mode 

— derivative with respect to time 
dt 

a angle of attack 

r qy cross-spectral coherence function for q(t) and y(t) 
p r effective radius of i& vibration mode 
o) r circular frequency of t*" mode 

° degrees 


Model Attitude Measurement Systems 


Inertial System 

The inertial AOA package used at LaRC [1] to sense the attitude of the model is shown 
installed in the nose of a test model in figure 1 . The AOA package uses a servo- 
accelerometer with its sensitive axis parallel with the longitudinal axis of the model. For 
quasi-static conditions, this sensor provides a model attitude measurement with respect to 
the local gravity field to an accuracy of ±0.01 ° over a range of ±20°. Most wind-tunnels use a 
cantilever arm (sting) to support the model. The model mounted at the end of the sting 
experiences dynamic oscillations due to unsteady flows that induce centrifugal accelerations 
on the inertial AOA package and result in angle of attack measurement errors. 

Physics of Problem. The physics of the problem are studied by considering the response of 
a single yaw mode as simple harmonic motion as depicted in figure 2. For the system shown 


3 



with natural frequency of oscillation, 
v r (f ) , can be written: 


co r , the AOA package displacement, y r (t ) , and velocity, 


y r {t) = Y r sm((o r t) 

v r (0 ■ = K cos ((o r t ) ; where V r = Y> f 


The corresponding tangential and normal acceleration components, a, and a n , are 

dv, 

a A f ) = ~^f m A sin((fy r f) ; where 4 = -Y r (o r 2 
a n(t) = ~ L ~ = ^-COS 2 (co r (t)) = ^-(1+cos(2o> r f)) 

Rewriting equation (4) in terms of peak acceleration gives 

A 2 

a n(t) = — ^ — (1 + cos( 2 to r fl) 

2 <*V Pr 


(3) 

(4) 

(5) 


The vibration-induced normal acceleration results in the AOA packaqe sensina a centrifunai 

acceleration coincident with its sensitive axis (i.e., the longitudinal axis of the model) The 
AOA package output prior to filtering becomes: or me model) . The 


V 4 

Anf(0 = fl'sina - ~-(1 + cos( 2 co r t)) + a x (f) 


( 6 ) 


The first term on the right-hand side of the equation is the gravitational acceleration due to 
die model attitude, a, relative to the local vertical. The second ter^sth^nS 
acceleration (from equation (4)) caused by the model yaw motion Accelerations resultino 

thTthlTt’erm 11 Tn l 0 /l 9 ! tu ^ nal model v ' bra,ions (typically high frequency) are represented by 
requereyT mean an 9 le ' ,he A0A signal is lowpass filtered (0.4 Hz cut-off 


. v r 2 

2 Pr 


( 7 ) 


It is important to note that the model vibration causes a bias error in the model attitude 
measurement that cannot be removed by filtering or data averaging. 
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Model pitch vibration causes a similar bias error term, where the tangential velocity is acting 
in the pitch plane. If the vibration response is composed of multiple yaw and pitch modes, 
the total bias error will be a linear summation of the error contributions for the m modes. 


a 

Afii » grsina - 2j~z~ 

r=r pr 


Or, in terms of the peak acceleration, from equation (5), 


( 8 ) 


Afif^gsma - 


r=1 2w r Pr 


(9) 


The above discussion is based on the case of continuous sinusoidal model motion. In the 
wind-tunnel, the data is non-stationary and random in nature. This results in a time varying 
bias error that is dependent on the number of modes participating and the amplitudes of 
motion for those modes (i.e., V r (t),A r (t),Afn(t)). 

Bias Error Correction. Several methods have been proposed to correct for this bias error 
using measured tangential accelerations at the AOA sensor location due to model yaw and 
pitch motion. One method is to measure the natural frequencies from the frequency spectra 
of the tangential accelerations and then determine the bias magnitudes by measuring the 
magnitude of the second harmonic components from the frequency spectrum of the unfiltered 
AOA signal (see equation (6)). This technique may be difficult to implement due to the 
participation of multiple modes and the required data accuracy to measure the small 
magnitudes at the second harmonic frequency. Young et. al. [7] determine the natural 
frequencies, w r , and corresponding peak acceleration magnitudes, A r , from the frequency 
spectra of the yaw and pitch acceleration measurements. The bias error is then calculated 
from the second term of equation (9). The vibration mode effective radius, p r , was 
estimated from sinusoidal input-output data taken during wind-off ground vibration tests and 
later confirmed with measured mode shapes. Both techniques use frequency domain signal 
processing which is suitable for the stationary data observed in the wind-off tests but is 
questionable when evaluating the non-stationary data observed during wind-tunnel testing. 

A third bias correction technique was developed and used at the National Aerospace 
Laboratory (NLR) in the Netherlands. This technique was developed for one vibration mode 
in each the yaw and pitch plane. Two additional accelerometers are used to measure the 
tangential accelerations due to the yaw and pitch motion of the model. The tangential 
accelerations are integrated to obtain velocity, squared, and divided by a scale factor to 
compensate for the effective radius of the vibration mode. This signal is then added to the 
unfiltered AOA output to cancel the second term of equation (6). The corrected signal is then 
lowpass filtered to obtain the corrected mean AOA measurement. The mode radius in the 
yaw and pitch plane is determined by tuning a potentiometer while manually exciting the 
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(3) Numerically integrate x(t), square, and denote by y(t); (4) Compute power spectra S qq 
and Syy, cross-spectrum, S qy , and cross-spectral coherence function r^; (5) Inspect T™, 
for spectral correlation within the AOA filter passband and identify the corresponding modal 
frequencies (denote as <o r ); (6) Estimate radius, p r , corresponding to each mode; 

(7) Band-pass filter lateral motion signal x(t) about modal frequency a> r and denote as x L (t ) ; 

(8) Low-pass filter dynamic AOA signal and denote by q^t ) ; (9) Numerically integrate and 
square x L (t) and denote by y^(f) ; (10) Correct inertial AOA output, qi_(t) -/£.(/)/ Pr - anc * 
lowpass filter the result; (1 1 ) Repeat for each mode. 


A sensitive correlation test between the time series q(t) and y(t) is provided by the cross 
spectral density coherence function, defined as follows: 


r , x SqyM 

qy( " )= 


( 10 ) 


where S q q(co) and Syy (to) are the power spectra of q(t) and y(t), respectively, and Sqy (<w) is 
the cross spectrum. Correlated spectral components common to both Sqq (to) and Syy (0)) 
appear in T q y(w) . Other spectral components common to S qq ((D) and Syyico) which are not 
phase coherent, i.e., unsynchronized, tend to be removed from Sqy (to) by averaging and 
cancelled by normalization and do not appear in T q y((o) . The coherence function and cross 
spectrum thus provide a means of detecting and quantifying AOA bias errors due to angular 
oscillation. 


Optical System 

A number of optical techniques have been investigated for measuring AOA [8-11], but none 
of these appear practical for both air and cryogenic (nitrogen) operations such as the NTF. 
Some of the problems with these and other optical techniques are that lasers are often used 
as light sources which require special considerations for packaging in a harsh wind-tunnel 
environment (especially a cryogenic high pressure environment such as the NTF). In 
addition, special sensors (some requiring wiring) may need to be placed in the model 
fuselage. Some of these devices may require windows in the fuselage which need special 
protection during model preparation. Another problem common to most of the optical 
techniques (including the video technique discussed here) is that fog in the test section may 
attenuate light to the point where the optical angle measuring system may not work at all. 

The onboard inertial sensor, however, is not affected by tunnel fog. Tunnel fog is not an 
operational problem in the NTF but may be encountered in other wind-tunnels. 

Recent success at NTF in measuring wing twist in both air and cryogenic runs led to the 
question of whether or not a similar technique could be used to measure model AOA at the 
facility. At NTF, wing twist is determined optically from photogrammetric measurements 
made on digital video images. The technique is based upon a single camera 
photogrammetric determination of two dimensional coordinates with a fixed (and known) third 
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dimensional coordinate. The wing twist is found from a conformal transformation between 
flow and no-flow 2-D coordinates in the plane of rotation. 

Since reduced data are not available in "real-time", the current wing twist technique was 
initially not deemed suitable as the primary means of determining AOA. However, it was 
thought that additional AOA data from such a technique (especially data at dynamic 
conditions) might aid in the analysis of the dynamic bias error problem associated with inertial 
sensors. For this reason, and to gain initial wind-tunnel experience with the technique for 
measuring AOA, data were taken for a number of runs over several days at both air and 
cryogenic conditions. Instead of using targets on the fuselage (in a manner similar to that 
used for wing twist), a fortuitous glint near the top of the fuselage caused by a specular 
reflection from one of the test section lights was used as a line source from which the slope 
and intercept could be computed and AOA determined. It was determined from this limited 
data set that if the onboard accelerometer was used for no-flow calibration at a particular 
tunnel temperature, then the AOA could be measured with flow to 0.02° root mean square 
(rms) in air mode and 0.05° rms in cryogenic mode when compared to the onboard 
accelerometer. Worst case disagreement with the onboard inertial sensor was 0.2 degrees. 

It was not determined during these tests how much of the disagreement was due to model 
dynamics causing bias errors in the inertial sensor measurements and how much of the 
disagreement was due to the video angle of attack measurement. Given the potential 
advantages of the video method, the results were considered encouraging enough to 
continue development and testing of the technique. 

Major advantages of the video angle measurement method compared to other optical 
techniques are that no laser is required since existing test section lights are used to illuminate 
the model, and no sensor is required in the model fuselage. In addition, for NTF, many of the 
operational problems associated with the packaging of instrumentation in a cryogenic high 
pressure environment, window frosting and equipment survivability have already been 
addressed during the development of the wing twist method. The difficulty of these 
operational problems should not be underestimated. Some of the solutions to these 
problems are not satisfactory, even after eight years of tunnel testing with video systems. 

One disadvantage of the image-based video technique is that lighting requirements for AOA 
may differ from the lighting requirements for general model surveillance. For very critical 
surveillance, the Test Engineer may find it necessary to re-adjust the test section lighting 
which might invalidate the video calibration or may even cause the video technique to be 
unable to reduce data at all. Another concern with the video method is the potential 
sensitivity of the technique to the location of the line to be measured on the image plane. 

This sensitivity to line image location may be difficult to completely compensate by 
calibration. Another disadvantage is that the technique does not function as an angle 
transducer to produce a voltage output proportional to angle as is common to most of the 
other optical techniques and inertial sensors. Wind-tunnel data acquisition systems are 
generally set up to handle such transducers in a standard manner. To measure angle with 
the video technique, computations on the gray scale of digital images are required. These 
computations, while amenable to automation for a fixed geometry and high contrast, make it 
difficult to automate the process to include variable geometry, especially if high contrast 
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cannot be maintained. In addition, camera alignment stability, practical in-place facility 
calibration, and various assumptions made for the single camera photogrammetric solution 
are all issues which have not been fully addressed. Further laboratory and tunnel tests are 
necessary to more completely characterize and develop the video technique before it is 
suitable for routine use. 

EXPERIMENTAL TEST SETUP AND PROCEDURES 

The test setup is shown in figure 5. Inertial and optical systems were used to measure the 
attitude of a cantilevered model system during dynamic tests under wind-off conditions. 

Instrumentation 

The model system was installed in a model assembly bay at the NTF as shown in figure 5. 
The mounting consisted of a "rigidly" supported cantilever (sting) that is positioned by a 
pitch-roll-translation mechanism. The model is attached to the sting through a six component 
strain gage balance. 

The model was instrumented with an inertial AOA package [1] maintained at a constant 
temperature of 160°F. The signal conditioner for the AOA package provides both an 
unfiltered (0-300 Hz bandwidth) signal and a filtered (0-0.4 Hz bandwidth) signal. Two 
miniature accelerometers were installed on the face of the AOA package to measure the 
accelerations tangent to the sensitive axis of the AOA sensor in the yaw and pitch planes of 
the model. In addition, four accelerometers were mounted external to the model fuselage to 
measure model yaw and pitch motion. 

An electrodynamic shaker was used to excite the model through a single point force linkage 
as shown in figure 6. The excitation was input in the pitch and yaw planes at the model 
fuselage hardpoints, approximately 36 inches forward of the balance moment center. Sine, 
modulated sine, and band-limited (5-60 Hz) random shaker inputs were used. 

Data Acquisition Systems 

Dynamic Signal Analyzer. A Hewlett Packard model 3566A dynamic signal analyzer (see 
figure 7) was used to provide the shaker stimulus and perform on-line time and frequency 
domain signal analysis. The 16 channel signal analyzer was used to monitor and record the 
shaker force input, model force balance outputs, AOA filtered (static) and unfiltered (dynamic) 
outputs, and model accelerations. This system was also used to monitor the model yaw and 
pitch moments which established the dynamic test conditions for acquiring model attitude 
measurements. 

Optical System. The optical technique used to determine model attitude for these tests is 
based upon the recording and analysis of digitized video images. Figures 8 and 9 show the 
optical system setup. A video signal from a standard RS-170 solid state camera with 
752 horizontal X 240 vertical pixels per field is routed to a frame grabber controlled by a 
486-66 MHz PC which records video fields every 1/60 second until the grabber memory 
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and the location of the cameras front perspective center was estimated from direct 
measurements. The camera was located 56 inches from the center of the fuselage. 
Normally the final calibration step is to run a no-flow pitch sweep over the range of angles of 
interest in which the onboard inertial sensor is used for angle calibration at each tunnel 
temperature. For the preparation bay tests, only differences from the no-load angle were of 
interest so that no angle calibrations were required. 


Data Acquisition System. The Data Acquisition System (DAS), shown in figure 10, is a 
12 channel digital data acquisition system with 1 6-bit resolution. For the tests, data were 
recorded at a rate of 800 samples per second per channel. Recorded channels included the 
unfiltered and filtered inertial AOA output, the auxiliary yaw and pitch accelerometer outputs, 
the six unfiltered balance components, the shaker excitation function, and an auxiliary optical 
AOA sensor output. 

Test Procedure 

The model was set at a prescribed angle of attack under static conditions. The model system 
natural frequencies were identified using sine sweep excitation of the model in the pitch and 
yaw planes. For each natural frequency of interest, a sinusoidal forced response test was 
conducted by controlling the shaker input amplitude to provide a defined peak-to-peak pitch 
or yaw moment on the model force balance. The control test variables were pitch moment 
for modes that had predominantly pitch motion, and yaw moment for modes that had 
predominantly yaw motion. Inertial and optical systems were used to measure the model 
attitude for a series of moment amplitude levels for sinusoidal excitation at a prescribed 
natural frequency of the model system. Sinusoidal forced response tests were conducted 
with the model set at pitch angles of 0°, 4.3° and 6°. 

In addition to the sinusoidal forced response tests, the model attitude was measured for 
modulated sine and random excitation tests. The modulated sine and random 
excitations/responses are more representative of the model dynamics observed in actual 
wind-tunnel tests. The majority of the modulated sine tests were conducted with a 0.25 Hz 
modulation of the first natural frequency in the pitch and yaw planes. The random excitation 
was only evaluated in the pitch plane for the model at a nominal angle of 0°. In each case, 
the inertial and optical systems were used to measure the model attitude for a series of 
moment amplitude levels. 


Discussion of Test Results 


Inertial System 

Modal Correction Technique. Attempts were made to excite three natural vibration modes in 
each of the yaw and pitch planes of the model system. The modes of interest were 
determined from previous modal surveys of the model system and are listed in Table 1 . The 
slight shift in mode frequencies indicated for mode 1 were attributed to differences in the 
assembly bay test setup (i.e., shaker attachment/alignment; boundary conditions to 
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anS ' e Change) - Measurements taken during previous wind-tunnel tests 

PMO H. mnl L T' ™ deS being excited wer « ^ approximately 8-" Hz and 
28-30 Hz [6] in the yaw plane. The radius for a vibration mode was determined usino a least 

nf thtT f,t |° f tH f m ° dal deformations on ^e model fuselage to project the rigid bodvmotion 
of the fuselage to a point of rotation (see figures 3 and 4). The vibration mode effective 

radius is then estimated as the distance from the mode rotation point to the inertial AOA 
sensor location in the mode, fuselage. These predictions indicated thS the Ss mav be 
positive or negative dependent on the vibration mode shape. For the case wera tha rad u .« . e 
posrtive, the point of rotation for the vibration mode is ahead of the AOA packaqe The 
significance of the sign of the radius is that the bias error effect mayTposte or negative 

upon Vlb u rat,0n mode being excited. This is demonstrated by the response for 
1 p y f. p ^ ne modes shown in figures 1 1 through 1 3. For the first two yaw modes (9 8 and 
1 . 9,5 H 2* the 'Indicated model angle change is negative when the model is beinq driven with 

2“ ati0n f t the m0de natural frepuenc y and then returns to itL 

when the shaker system is shutoff. The third yaw mode (29.8 Hz ), which has a neaative 

radius value, shows an indicated positive angle change when the model is beinq driven with 
whPnThp rj tat ' 0r \ at the m ° de natUral frepuenc V aad then returns to te 

4.p U g^v^ 

structure resulting in high dynamic yaw moments with energy in the 28-30 Hz band This 
points out the need to do modal testing with the model installed in the uinnel 

Table 1 

Model System Vibration Mode Characteristics 


Mode No. 

Natural 

Frequency 

— m , 

Effective 

Radius 1 

(inches) 

Mode Description 

1 

8.8 to 9.8 2 

31.0 

Stinq First Bendina in Yaw 

2 

9.2 

30.2 

Stinq First Bendina in Pitrh 

3 

19.5 

0.18 

Model Yawing on Balance with Sting 
First Bendinq 

4 

r 20.3 

-1.08 

Model Pitching on Balance with Sting 
First Bendinn 

5 

29.8 

-7.14 

Model Yawing on Balance with Sting 
Second Bendinq 

6 

Note- I1\ Tb 

34.2 

-7.65 

x xi 

Model Pitching on Balance with Sting 
Second Bendina 


. . ; v H 11 Ul ,uiai|on ror a mode to the AOA sensor inratinn 
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The results of sinusoidal excitation tests (model at nominal angle of 0°) for the first mode in 
the yaw and pitch plane is shown in figures 1 4 and 15. For a set excitation level, time 
domain data were acquired and stored using the Hewlett Packard signal analyzer. These 
data were transferred to a personal computer where a software routine, written as an M-file in 
the MATLAB [14] language, was used to estimate the bias error in the inertial device. The 
yaw and pitch plane accelerations were bandpass filtered, integrated to obtain velocity, 
squared, divided by the mode effective radii, and scaled to obtain the equivalent angle 
change in degrees. The estimated bias error is in good agreement with the indicated mean 
angle change measured with the onboard inertial AOA sensor. If the estimated bias error is 
subtracted from the indicated model angle, the error would be reduced from a maximum of 
-0.146° to -0.009° for the first mode in the yaw plane (y) and from -0.175° to -0.006° in the 
pitch plane (z). These corrections are within the AOA accuracy requirement of 0.01 °. Similar 
results were obtained for the sinusoidal input tests with the model at nominal angles of 4.3° 
and 6°. 

In addition to the sinusoidal tests, the bias error was examined for modulated sine and 
random inputs. Figure 1 6 shows the estimated and measured bias error as a function of 
time for a 9.2 Hz pitch excitation with a 0.25 Hz modulation. Excellent agreement is obtained 
with the difference between the estimated and measured bias being less than 0.005°. 
Modulated sine tests were conducted for the first mode in each the y and z axes at several 
excitation amplitude levels and consistent results were obtained between the measured and 
predicted bias errors for all cases. In the pitch plane, two levels of random excitation were 
input to the model. An eight second record of the inertial AOA sensor response for the 
highest level random excitation is shown in figure 17. The random response measured by 
the vertical (pitch plane) accelerometer on the face of the AOA package was composed of 
primarily 9.2 Hz response. The bias estimate based on only the 9.2 Hz mode contribution is 
also shown in figure 17. Again, the estimated and measured bias are in very good 
agreement. 

Time and Frequency Domain Analysis. The utility of spectral techniques in identifying 
dynamic AOA data records corrupted by angular motion is illustrated in the following data set 
obtained during yaw excitation. The model is excited in yaw at 5300 inch-pound amplitude at 
9.75 Hz, modulated at 0.3 Hz. Figure 18 illustrates a 70 second record of the filtered AOA 
time series and the corresponding integral-squared yawing moment. Figure 19 depicts their 
corresponding power spectra. The apparent good agreement at 0.3 Hz seen both in the 
time series and in the power spectra is corroborated by the well defined 0.3 Hz peak in the 
cross spectral coherence shown in figure 20, and by the absence of significant energy below 
0.3 Hz in both power spectra. AOA bias correction to better than 0.01 ° was obtained in the 
first third of the record. 

For comparison, the integral-squared y-axis accelerometer time series and its power 
spectrum appear in figures 21 and 22. Note that the time series agreement with AOA is poor 
in the latter third of the record, and note the presence of significant power spectral energy 
from 0 to 0.3 Hz in the integral-squared y-axis accelerometer record, which is absent in the 
AOA power spectrum. The cross spectral coherence seen in figure 23 shows no correlation 
below 0.3 Hz and less pronounced correlation at 0.3 Hz than for yaw. Since the low 


13 



SSSiSSS?* r 
^issssssr ■ - — — * ttzxzssszfx 

procedure was applied tothe^odulMedtin^ 

z^iK^ U terOT^ter' n Sete!rteTresults^re S ^ttetedtnT^le2 m h nl "“T' ,0r ° e ' 
aSTndtrde“rn1f'r f rS urK^ec^°AOA -an 

regression, the mean of ,he corrected fn"^ 

Table 2 

Estimated Bias Error with Pitch Excitation 
(Unexcited Data Zero, AOA = -0.290°) 
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Single mode bias correction by multiple linear regression using the previously described 
frequency and time domain procedure was also applied to the modulated sine data sets for 
yaw excitation. Selected results are tabulated in Table 3 showing the time and duration of 
the corrected time series, the source of the integral-squared correction, the mean and 
standard deviation of the filtered uncorrected AOA, the bias error estimated by the 
regression, the mean of the corrected AOA, and the standard error of the regression. 

Table 3 

Estimated Bias Error with Yaw Excitation 
(Unexcited Data Zero, AOA = -0.283°) 


Record 

Segment 

(Sec) 

Correction 

Source 

Mean 

Filtered 

AOA 

(Deg) 

Standard 

Deviation 

Filtered 

AOA 

(Deg) 

Bias 

Estimate 

(Deg) 

Mean 
Corrected 
AOA (Deg) 

Standard 

Error 

(Deg) 

10-30 

Yaw 

-0.433 

0.078 

-0.148 

-0.285 

0.007 

! 40-60 

Yaw 

-0.440 

0.077 

-0.161 

-0.280 

0.032 

i 10-30 

Side 

-0.433 

0.078 

-0.178 

-0.254 

0.050 

10-30 

Y Accel. 

-0.433 

0.078 

-0.149 

-0.284 

0.021 


Figure 26 illustrates the residuals of the AOA correction using integral-squared yaw 
computed over 10 to 30 seconds. Note that the AOA error after correction is reduced by 
91 percent. In Table 3, the correction uncertainty increases tenfold in the 40 to 60 second 
interval because of a significant low frequency random disturbance in the integral-squared 
yaw which is absent in the AOA time series. 

Optical System 

The purpose of these tests was to determine the effect of dynamics on the video angle 
measurement method and compare with results for an onboard inertial sensor. For these 
dynamic tests the change in model AOA was of primary interest rather than the measurement 
of the AOA itself. In order to measure the error in indicated angle due to dynamics, it was 
necessary to determine the basic instrumentation repeatability. This repeatability was 
measured as the standard deviation of data recorded over time with the pitch angle fixed and 
could be classified as either short term (1 sec) or long term (several days). 

The short-term repeatability for the video technique over 1 second (60 fields) was computed 
as the mean of the standard deviation at each of 36 “shaker off data points and found to be 
0.0018 ° with a maximum value of 0.0034 °. The mean standard deviation of the y-intercept 
for the “shaker off points over 1 second was 0.0002 inch with a maximum value of 
0.0006 inch. 

The shaker was attached to the fuselage for measurements in yaw for the first 42 data points 
of the test, which occurred over a 48-hour period. Since the model pitch angle was assumed 
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unchanged during this period, the variability in the mean pitch angle measured for "shaker 
Ch “f s ,he lon 9- ,erm repeatability of the angle 4asunSeL ?he 
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Table 4 

Summary of AOA Error for Video and Inertial AOA Measurements 


Technique 

Mean 

Standard 

Deviation 

Maximum 

No. of Data 
Points 

Excitation 

Video 

0.000° 

0.002° 

0.006° 

n 49 

Yaw 

Inertial 

-0.050° 

0.057° 

! -0.217° 

49 

Yaw 

Video 

0.002° 

! 0.004° 

0.010° 

42 1 

Pitch 

Inertial 

-0.033° 

0.046° 

-0.174° 

42 

Pitch 


large, generally negative, bias error that can be as larae asn? Henroao q e er ® uffers a 

«.e AOA error in yaw and pitch for the vid^and ^ 

to-peak moment in figure 27, AOA errors for the various runs are plo^S 

ynamrc yaw tests were conducted at several model pitch angles, whereas, the dynamic 
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pitch tests were limited to pitch angles near 0 degrees because of shaker mounting 
constraints. 

Selected temporal variations over one second in the video measured angle during shaker 
excitation in yaw at 9.8 Hz are plotted in figure 31 . The video technique recorded changes in 
pitch during yaw excitation of up to 0.04° at the largest yaw moment of almost 4000 in-lb 
peak-to-peak. The corresponding y-intercept plots are presented in figure 32 where 
displacements in the pitch plane of up to 0.01 inch are noted for the largest yaw moment. 
Example temporal angle variations over 4.6 sec during shaker excitation in pitch are plotted 
in figures 33 and 34 for 9.2 Hz sine excitation modulated at 0.25 Hz, and band-limited 
(5-60 Hz) random excitation. Pitch variations of almost 0.3 ° were noted during excitation of 
the 9.2 Hz pitch mode. The corresponding y-intercept plots are presented in figures 35 and 
36 which show displacements in the pitch plane of up to 0.1 inch for the 9.2 Hz mode. 

The standard deviation of the video measurements over 1 second versus yaw and pitch 
moment are presented in figure 37. The maximum standard deviation over 1 second for yaw 
moment is generally less that 0.02°, even for the largest moments, indicating very little 
dynamic pitching during yaw dynamics. The maximum value for the largest pitch moment is 
10 times larger due to the dynamic changes in pitch angle during dynamic pitch. 

Concluding Remarks 

Servo-accelerometers are excellent devices and widely used for measuring angle of attack 
for wind-tunnel models. However, when these inertial devices are used in the presence of 
significant model system vibrations, bias errors are introduced due to centrifugal 
accelerations associated with model system oscillations. This bias error can be difficult to 
estimate and correct. An experimental study of inertial and optical model attitude 
measurement systems response in a wind-off dynamic environment is presented in this 
paper. This study is a benchmark test for developing new model attitude measurement 
systems as well as studying techniques for correcting existing inertial sensors in order to 
obtain high accuracy angle of attack measurements in a hostile dynamic environment. The 
results of this study will be used for evaluating model attitude measurements using the 
inertial and optical systems during wind-tunnel tests of the model system testbed. Further 
work needs to be done to extend the AOA correction algorithms to multiple mode 
excitation/response and validate the algorithms during actual wind-tunnel tests. 

The modal correction technique developed and used for calculating the AOA bias error 
worked very well for the wind-off dynamic environment of sinusoidal, modulated sinusoidal, 
and random forced response testing. It was found that some of the model system vibration 
modes excited to high levels during wind-tunnel tests in the NTF could not be excited very 
well in the NTF model assembly bay. This is attributed to boundary conditions at the 
assembly bay model support which are quite different from the NTF model support system. 
This points out the need for doing modal testing on the model system after installation in the 
tunnel to assure the most accurate measurement of the vibration mode radii and for the 
evaluation of the AOA response to model system vibrations. 
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Figure 1.- Schematic of wind tunnel model system. 


Figure 3.- Sting bending in yaw plane, 
9.0 Hz vibration mode. 
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Figure 4.- Model yawing on balance, 
29.8 Hz vibration mode. 
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Figure 2.- EWectof vibration on Inertial modal 
attitude measurement. 
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Figure 5. Test setup in model assembly bay. 
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Figure 10. Data acquisition system. 
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Figure 1 1 . Inertial AOA measurement, yaw 

acceleration, and yaw moment versus 
time for 9.0 Hz sinusoidal input in yaw 
plane. 
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igure 12. Inertial AOA measurement, yaw 

acceleration, and yaw moment versus 
time for 19.5 Hz sinusoidal input in yaw 
plane. 
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Figure 13. Inertial AOA measurement, yaw 

acceleration, and yaw moment versus 
time for 29.8 Hz sinusoidal input in yaw 
plane. 



Figure 14. Measured and estimated AOA bias 

error versus peak-to-peak yaw moment 
for 9.0 Hz yaw mode. 



Figure 15. Measured and estimated AOA bias 
error versus peak-to-peak pitch 
moment for 9.2 Hz pitch mode. 



Figure 16. Measured and estimated AOA bias 

error versus time for 9.2 Hz sinusoidal 
excitation in pitch with 0.25 Hz 
modulation, maximum peak-to-peak 
pitch moment of 4500 inch-pounds. 



Figure 17. Measured and estimated AOA bias 

error versus time for random excitation 
in pitch. 
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7-Frm Avg Spectrum of AOA Dynamic 



7-Frm Avg Speclrum of lnl(Yaw Momanl )*2 



Figure 19. Power spectra of AOA output and integral-squared yawing moment for yaw excitation 
of 5300 inch-pounds at 9.75 Hz with 0.3 Hz modulation. 

Cross-Spectraf Coherence 



Figure 20. Cross-spectral coherence between AOA output and integral-squared yawing 

moment for yaw excitation of 5300 inch-pounds at 9.75 Hz with 0.3 Hz modulation. 
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Yaw points only Pitch points only 
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27. Indicated AOA bias errror versus peak-to-peak moment for video method (denoted with 
circles) and onboard inertial AOA device (denoted with X’s) for all yaw and pitch tests. 
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Figure 28. Comparison of AOA bias error versus peak-to-peak yaw moment for video method (denoted 
with circles) and onboard inertial AOA device (denoted with X’s) with the excitation frequency, 
modulation frequency (in parantheses where appropriate), and nominal angle identified at the 
top of each plot. 
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Figure 29. Comparison of AOA bias error versus peak-to-peak yaw moment for video method (denoted 

n? > and on .b oar d mert'ai AOA device (denoted with X’s) with the excitation frequency 
and nominal angle identified at the top of each plot. q Cy 
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Figure 30. 
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Figure 31. Change in pitch angle versus time measured with the video method for 9.75 Hz yaw excitation 
at yaw moments of 3983 in-lbs and zero. 




Figure 32. Corresponding change in Y-intercept versus time as measured with the video method for the 
same measurement points as Figure 31 . 

.4 481 Pitch 9.2 HZ (.25 Hz nod) 4509 ln-lb 
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Tint, sec 

Figure 33. Change in pitch angle versus time 

measured woith the video method for 
9.2 Hz pitch excitation with 0.25 Hz 
modulation and a maximum peak-to- 
peak pitch moment of 4509 in-lbs. 


Tina. aac 

Figure 34. Change in pitch angle versus time 

measured with the video method for 5- 
60 Hz random excitation in pitch plane. 
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Figure 35. Corresponding change in Y-intercept 

versus time as measured with the video 
method for the same measurement 
point as Figure 33. 
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Figure 37. Video method pitch angle standard deviation 
pitch moment of the model. 
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Figure 36. Corresponding change in Y-intercept 

versus time as measured with the video 
method for the same measurement 
point as Figure 34. 



over one second versus peak-to-peak yaw and 


38 





Form Approved 

i OMBNo. 0704-0188 

Public reporting burden for this objection of information is estimated to average 1 hour rw ZZ33HI ; oj a 

r^- r;: va rr • and ,o — - 


2. REPORT DATE 

February 1995 


1 4. TITLE AND SUBTITLE 

Dynamic Response Tests of Inertial and Optical Wind-Tunnel Model 
Attitude Measurement Devices 


3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 

5. MINDING NUMBERS 


1 6. AUTHOR(S) 

R. D. Buehrle, C. P. Young, Jr., A. W. Burner, J. S. Tripp, P. Tchena 
T. D. Finley, and T. G. Popernack, Jr. 

1 7. PERFORMING ORGANIZATION NAME(S) AND ADDRESSES) 

NASA Langley Research Center 
Hampton, VA 23681-001 


1 9. SPONSORING / MONITORING AGENCY NAME(S) AND AODRESS(ES) 

National Aeronautics and Space Administration 
Washgton, DC 20546-0001 


WU 505-59-54-01 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


10. SPONSORING / MONITORING 
AGENCY REPORT NUMBER 

NASA TM-109182 


|l1. SUPPLEMENTARY NOTES 

TT”' If. 6 "®' Finley ’ Po ' ,ernack: Langley Research Center. Hampton. VA 
Young. North Caroline State University, Raleigh, N.C. H 


|12a. DISTRIBUTION / AVAILABILITY STATEMENT 

Unclassified - Unlimited 


Subject Category 09 


12b. DISTRIBUTION CODE 


|13. ABSTRACT (Maximum 200 words) 

!S?f are presen,ed for an experimental study of the response of inertial and optical wind-tunnel model 
attrtude measurement systems in a wind-off simulated dynamic environment This studv is Dart of aTomo,™ 

SS2S.S 1 ' Cemerto^h^accuracy.^^SefruSrrin, 

systems that can be used in a dynamic wind-tunnel environment. This activity was Dromnted hv the inoniai 

i£ 

degree of accuracy forTheTbS The oplical 

system results were uncorrupted by model vibration in the laboratory setup The ° pt,Cal 


|14. SUBJECT TERMS 

Model Attitude Measurement; Wind Tunnel; Model Vibration; and Bias Error 


117. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


NSN 7540-01-280-5500 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 


IS. NUMBER OF PAGES 
43 


16. PRICE CODE 
A03 

20. LIMITATION OF ABSTRACT - 


Standard Form 298 (Rav. 2-89) 

Praannkwl Ki/ AMcm Piw Tnn 4 - ' 


